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In this talk, we briefly review the current understanding of the behav- 
ior of the QCD equation of state throughout the phase diagram. Special 
emphasis is given to regions of phenomenological interest, and a number 
of important open questions as well as directions of ongoing research are 
pointed out. These include in particular the region of low temperatures 
and (moderately) high densities, where at the moment we have extremely 
few first principles tools available. 

PACS numbers: ll.lO.Wx, 12.38.Mh, 21.65.Mn, 21.65. Qr. 

1. Introduction 

Apart from being of obvious theoretical interest, the behavior of the 
Equation of State (EoS) of QCD as a function of temperature T and baryon 
chemical potential is of considerable phenomenological importance in 
a variety of physical systems. Through Einstein's equations, the high- 
temperature and low-density EoS governed the expansion and cooling rates 
of the universe when its age was of the order of a few tens of microseconds, 
and through the expansion rate it thus had an effect on e.g. the decoupling 
times of various dark matter candidates. In heavy ion physics, the success 
of the hydrodynamic description of the collision relies on the availability of 
accurate information on the energy density and pressure of the underlying 
theory. And finally, the properties of compact stars — in fact, their entire 
composition — are extremely sensitive to the EoS of cold and dense nuclear 
and quark matter. 

In this talk, we will review what is currently known about the EoS of 
QCD in various parts of the fiB-T phase diagram, covering both the confined 
and deconfined regimes and extending from zero density all the way to zero 
temperature. Conceptually, the setup of the problem is extremely simple: 
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One takes the Euclidean functional integral corresponding to the partition 
function (or grand potential) of the theory, 

n{T,{nf},{mf}) = -riog j Vm'I^A^e--^o'^^^'^''''^QCB^ (1) 
-Cqcd = \F;,F;, + i^f{-i^D^ + mf - iif^o)i^f, (2) 

evaluates it one way or another, and subsequently uses standard relations to 
obtain predictions for a number of equilibrium thermodynamical quantities 
ranging from the pressure to the quark number, entropy and energy densities 



pv = -n, (3) 

= -df.f^^ (4) 

sV = -5tO, (5) 

e = -p + Ts + fifUf. (6) 



The real challenge of course becomes, how to choose the optimal method 
for doing this in different parts of the phase diagram, and this question is 
what we will be mostly concerned with. There is no one single method to 
cover the entire ^_b-T plane, but as we will argue below, combining results 
from various existing techniques provides us with a quantitative handle on 
the behavior of the EoS almost everywhere. 

To make the following discussion more transparent and physically in- 
tuitive, we choose to parametrize the temperature and baryon chemical 
potential by a radial and angular variable. These are defined by 

T = pcosO, jiB = \/67rpsin^, (8) 

where p is roughly proportional to the fourth root of the energy density of 
the free theory, while the angle measures the deviation from the temper- 
ature (zero density) axis. 

The roles of the variables are clear: Due to the running of the , the 
energy scale p measures how strongly coupled the system is. Roughly speak- 
ing, at p < 100 MeV we are deep in the confined phase of the theory, at 
100 McV p 1 GeV the system resides in an intermediate region charac- 
terized by a deconfinement phase transition (or crossover) and complicated 
dynamics, and at p > 1 GeV, we enter a regime where a description in terms 
of weakly interacting, deconfined quasiparticles becomes feasible. The angle 
d on the other hand separates two phenomenologically interesting regions: 
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The values <1 are relevant for the physics of the early universe and heavy 
ion collisions, while in the matter found inside neutron stars, we have to a 
very good accuracy 9 = tx 12. Thus, while it may be theoretically interesting 
to consider the EoS for other values of Q as well, most of our attention will 
in what follows be divided between these two cases. 

2. Small Q: Prom Hadron Gas to Quark-Gluon Plasma 

For baryon chemical potentials of at most the order of the temperature, 
the two limits of low and high energy densities, respectively p < 100 MeV 
and /c > 1 GeV, are conceptually relatively simple. At low energies, we 
are dealing with a dilute gas of hadrons that may be described in terms of 
hadron resonance gas (HRG) models, while at high p, the system comprises 
a (relatively) weakly coupled quark-gluon plasma that is amenable to treat- 
ment via resummed perturbation theory. In the regime between these two 
limits, one may attempt to gain qualitative understanding of the dynamics 
of the system via various types of effective theories and models, but the 
ultimate tool for quantitative information is clearly lattice QCD. 

The success of the hadron resonance gas models, originally proposed in 
Ref. [1], is based on confinement. At sufficiently low energy density, the sys- 
tem comprises a dilute gas of hadrons, in which the strong interactions are 
very efficiently confined inside these particles, making the system in practice 
a free Stefan-Boltzmann gas. Raising the value of p, and thus introducing 
more degrees of freedom, one must begin to take into account both the in- 
teractions between the hadrons and the fraction of spatial volume occupied 
by them, usually according to the Van der Waals prescription. Eventually, 
these effects lead to sizable uncertainties in the results. Nevertheless, cur- 
rent calculations, which typically include O(IOO) hadrons up to masses of 
2-3 GeV, have been extended to close to the phase transition region, and 
thus the results now have a finite overlap region with lattice simulations. 
For details of these calculations as well as up-to-date results, see e.g. Ref. [2] 
and references therein. 

For p between roughly 100 MeV and several times the critical tempera- 
ture Tc of the deconfinement transition, the method of choice for quantita- 
tively reliable results is lattice QCD. For a long time, the main challenge in 
lattice simulations was to reach physical quark masses and a wide enough 
temperature range, but recent years have witnessed enormous progress in 
both of these fronts [SlH]. As demonstrated e.g. in Fig. 2 of Ref. [5], there 
is now an impressive agreement between current lattice data and the pre- 
dictions of hadron resonance gas models at temperatures slightly below T^. 
In addition, the discrepancies between the predictions of the different lat- 
tice groups regarding the EoS and the critical temperatures of the chiral 
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and deconfinement transitions have recently significantly decreased. At the 
moment, the main limitation of the lattice approach to QCD thermodynam- 
ics is clearly the sign problem, prohibiting simulations at baryon chemical 
potentials much larger than the temperature. For recent progress in this 
direction, see e.g. Ref. [6j and references therein. 

Even higher temperatures and energy densities, p ^ 1 GeV, witness 
lattice simulations becoming more and more expensiv^i] while perturbation 
theory starts to work better and better. While the pure weak coupling ex- 
pansion of the QCD pressure, at the moment worked out up to the partial 
term [HI [9l [lOl [11], shows relatively poor convergence, the situation can be 
greatly improved using resummed perturbation theory. The most natural 
platform for this is offered by dimensionally reduced effective theories, in 
addition to which very promising results have recently been obtained using 
Hard Thermal Loop resummed perturbation theory |121ll3j. The perturba- 
tive approach has in addition the virtue of being easily generalizable to finite 
chemical potentials [14], as even dimensional reduction has been shown to 
lead to quantitatively reliable results as long as f^B ^ T/g [ T^ . 

3. 6» ?i; n/2: Cold Nuclear/Quark Matter 

The region of the QCD phase diagram characterized by finite baryon 
density and (negligibly) small temperatures is in many ways more problem- 
atic than that of small 9. For one thing, lattice QCD is no longer available 
as a computational method, and one must thus base one's approach on a 
combination of weak coupling techniques and model calculations. In addi- 
tion, it appears that the exact phase structure of the theory has a highly 
non-trivial form at small temperatures, and it is only at asymptotically large 
chemical potentials that the true nature of the deconfined phase is known. 
There, the pairing instability of the perturbative ground state drives quarks 
on the Fermi surface to pair in a way that ties their color and flavor indices 
together, resulting in the so-called Color Flavor Locked (CFL) superconduc- 
tor. The questions of how far down the fiB axis this phase extends to and 
which of the several proposed candidate phases ultimately replaces it, are, 
however, still open. Luckily, not all equilibrium thermodynamic quantities, 
in particular the EoS, are extremely sensitive to these details. 

Starting again from small energy densities and the nuclear matter phase, 
we enter a corner of the phase diagram that is under relatively good quan- 
titative control. This is not due to first principle calculations, however, but 
rather to the abundance of experimental data on low energy nucleon-nucleon 
scattering, enabling an accurate modeling of the dilute nuclear matter EoS. 
In close analogy to the HRG approach, the quantitative accuracy of these 



^ See, however, Ref. [7] for lattice results extended to high temperatures. 



excqcdvuorinen printed on January 20, 2013 



5 



models begins to suffer as one increases p and approaches the deconfinement 
transition region. In the case of cold and dense nuclear matter, the primary 
uncertainties in the EoS are due to the unknown composition of the matter 
(the conjectured presence of e.g. hyperons or kaon condensation), and mi- 
nor ones to the details of the model calculations such as the exact form of 
the respective variational ansatz or the effect of neglecting the simultaneous 
interactions of more than two nuclei. For further details of these computa- 
tions as well as for references regarding the nuclear matter EoS, the reader 
is asked to consult Refs. [i mfTTlfTH] . 

In the case of asymptotically large energy densities, p » 1 GeV, the 
physical ground state of QCD is a CFL superconductor. It is characterized 
by a (perturbatively) exponentially small energy gap A ~ e~*/^, and it is 
thus natural that the effects of the gap can be formally neglected in a weak 
coupling expansion of the EoS. In any practical applications, e.g. when dis- 
cussing compact stars, this is, however, certainly not good approximation, 
and one should in addition pay close attention to quark mass threshold 
crossings. Both of these effects have been included in the recent three- 
loop perturbative computation of Ref. |19j . The results of this calculation 
demonstrate that while the zero-temperature pressure does exhibit some- 
what better convergence properties than the high-T one, the uncertainties 
in the results become sizable when one approaches densities that might be 
realized in nature. 

It should not come as a surprise that in the ~ 7r/2 case, the most 
problematic region is by far that of intermediate densities, 100 MeV ^ p ^ 
1 GeV. In this regime, neither the nuclear model approach nor perturbation 
theory is applicable, but obtaining accurate information on the behavior 
of the EoS is nevertheless crucial in order to understand the structure of 
compact stars. Two very different lines of attack have been suggested to 
overcome these difficulties. One may either trust one of the several models 
that have been built to describe the relevant physics in this regime (usually 
following the guidance of symmetry principles) [20], or interpolate between 
the trusted EoSs of dilute nuclear and dense quark matter [19]. Recent 
advances in compact star observations have in addition made it possible to 
directly rule out sizable regions of their mass-radius plane. One day, this 
may even lead to the solution of the inverse problem, i.e. predicting the 
EoS of cold nuclear and quark matter solely from experimental data (see 
e.g. Refs. [211 [22] for a discussion of this). 
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